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Abstract — The relation between parity-check matrices of quasi- 
cyclic (QC) low-density parity-check (LDPC) codes and biadja- 
cency matrices of bipartite graphs supports searching for power- 
ful LDPC block codes. Using the principle of tailbiting, compact 
representations of bipartite graphs based on convolutional codes 
can be found. 

Bounds on the girth and the minimum distance of LDPC 
block codes constructed in such a way are discussed. Algorithms 
for searching iteratively for LDPC block codes with large girth 
and for determining their minimum distance are presented. 
Constructions based on all-ones matrices, Steiner Triple Systems, 
and QC block codes are introduced. Finally, new QC regular 
LDPC block codes with girth up to 24 are given. 

Index Terms — LDPC code, convolutional code, Tanner graph, 
biadjacency matrix, tailbiting, girth, minimum distance 



I. Introduction 

Low-density parity-check (LDPC) codes, invented by Gal- 
lager [1] in the 1960s, constitute a hot research topic since they 
are a main competitor to turbo codes (2|-|5l. Recently, a con- 
nection between LDPC codes and codes based on graphs was 
shown (see, for example, |6|-[10|), which opens new perspec- 
tives in searching for powerful LDPC codes. Moreover, coding 
theory methods can be applied in describing and searching for 
graphs better than previously known. For example, in ifTTI . 
1 12 1 compact representations based on convolutional LDPC 
codes for famous bipartite graphs such as Heawood's, Tutte's, 
and Balaban's graphs |fl3ll are presented. 

Typically, LDPC codes have a minimum distance which 
is less than that of the best known linear codes, but due to 
their structure they are suitable for low-complexity iterative 
decoding, like for example the believe-propagation algorithm. 
An important parameter determining the efficiency of iterative 
decoding algorithms for LDPC codes is the girth, which 
determines the number of independent iterations [1| and is 
a parameter of the underlying graph. The minimum distance 
seems not to play an important role within iterative decoding 
algorithms, since the error-correcting capabilities of such a 
suboptimal procedure are often less than those guaranteed by 
the minimum distance. In fact, it was shown in |14| that 
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the performance of LDPC codes in the high signal-to-noise 
(SNR) region is predominantly dictated by the structure of 
the smallest absorbing sets. However, as the size of these 
absorbing sets is upper-bounded by the minimum distance, 
LDPC codes with large minimum distance are of particular 
interest. 

LDPC codes can be characterized as either random/pseudo- 
random or nonrandom, where nonrandom codes can be subdi- 
vided into regular or irregular Q, O, flj], HSJ-EU, while 
random/pseudo-random codes are always irregular (29], 1 30 1 . 
A (J, if) -regular (nonrandom) LDPC code is determined by a 
parity-check matrix with exactly J ones in each column and 
exactly K ones in each row. 

The class of quasi-cyclic (QC) (J, if)-regular LDPC codes 
is a subclass of regular LDPC codes with low encoding com- 
plexity. Such codes are most suitable for algebraic design and 
are commonly constructed based on combinatorial approaches 
using either finite geometries lfl5l or Steiner Triple Systems 
lfl6l . II 171 . having girth g > 6. Amongst other algebraic 
constructions leading to QC LDPC codes with larger girth 
we would like to mention |fl9l , where a class of QC LDPC 
codes of rate R = 2/5 with girth up to 12 based on subgroups 
of the multiplicative group of the finite field F p was obtained. 
The same method was used for convolutional codes in Ell . 

Although QC LDPC codes are not asymptotically optimal 
they can outperform random or pseudorandom LDPC codes 
(from asymptotically optimal ensembles) for short or moderate 
block lengths [ 18 1. This motivates searching for good short QC 
LDPC codes. 

The problem of finding QC LDPC codes with large girth 
and large minimum distance for a wide range of code rates 
was considered in several papers, for example, Q, 0, [18|, 
[21 1, [23 1, 1 24 1 . Codes with girth at most 12 are constructed in 
COD, ED, ED, ED, while O gives examples of rather short 
codes with girth 14. Codes with girth up to 18 with J > 3 
are presented in ll24l and it is shown that QC LDPC codes 
with girth > 14 and block length between 34,000 and 92,000 
outperform random codes of the same block length and rate. 

Most of the papers devoted to constructing nonrandom 
LDPC codes with large girth combine some algebraic tech- 
niques and computer search. Commonly these procedures start 
by choosing a proper base matrix (also called weight or degree 
matrix) or the corresponding base graph (also called seed 
graph ||20l or protograph QTl ). The references [18|, [20| are 
focused on all-ones base matrices, while in [7|, [24] base 
matrices are constructed from Steiner Triple Systems and 
integer lattices. In both cases, a system of inequalities with 
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integer coefficients describing all cycles of a given length 
is obtained and suitable labels or degrees are derived. For 
example, if we replace all nonzero entries in the base matrix 
by permutation matrices ll27l . IT3T1 . circulant matrices [8|, 1 18 1, 
ED, ED* ED, or sums of circulant matrices 11251 . we obtain 
the corresponding QC LDPC block codes. On the other hand, 
if we replace all nonzero entries in the base matrix by either 
monomials or binomials, we obtain the corresponding (parent) 
LDPC convolutional codes EYl. ll25l. Il28l. 

Notice that both constructing the inequalities and the label- 
ing require significant computational efforts. Some methods 
directed towards reducing the computational complexity of 
these steps can be found in l23l , |26|. 

Parameters of the so far shortest QC LDPC block codes 
with J = 3 and girth 6, 8, and 10 found via computer search 
are presented in [26|, improving previous results from (18). 

As mentioned earlier, it is important that the constructed QC 
LDPC block codes have large minimum distance for achieving 
a suitable upper-bound on their error-correcting performance 
at high SNR. It is proved in [2| that the minimum distance 
of QC LDPC codes whose base parity-check matrices are 
JxK all-one matrices is upper-bounded by (J+l)!. However, 
considering base matrices with zeros leads to QC LDPC codes 
with larger minimum distance. For example, in [25] it is 
shown that replacing all nonzero entries in the base matrix 
by sums of circulants and all zero entries by all-zero matrices, 
increases the minimum distance of the resulting code while 
preserving its regularity. For LDPC convolutional codes this 
approach implies that a parity-check matrix contains binomials 
instead of monomials. The corresponding upper-bound on the 
minimum distance of such LDPC codes is presented in l25ll . 
A particular case of this upper-bound, valid only for codes 
with zeros and monomials is derived in ifTTI . 

In Section HU we introduce notations for generator and 
parity-check matrices of convolutional codes and for their 
corresponding tailbiting block codes. Section [III] focuses on 
bipartite graphs, biadjacency matrices, and their relation to 
parity-check matrices of LDPC block codes. Our constructions 
of base and voltage matrices, used when searching for LDPC 
block codes with large girth, are introduced in Section [IV] 
Bounds on the girth and the minimum distance for QC (J, ir- 
regular LDPC block codes are discussed in Section [V] New 
search algorithms for QC LDPC block codes constructed 
from all-one matrices, Steiner Triple Systems, and QC regular 
matrices are presented in Section [VI] Moreover, depending 
on the desired girth, algorithms of different complexity for 
constructing the set of inequalities and searching for suitable 
labelings are described. A new algorithm for computing the 
minimum distance of QC (J, if) -regular LDPC codes is 
described in Section IVIII and used to compute the minimum 
distance of our newly found codes. Moreover, we determined 
the hitherto unknown minimum distance for some of the 
shortest known LDPC codes given in [24 1. In Section [Villi we 
present new examples of ( J, K )-regular QC LDPC codes in 
the form of tailbiting LDPC codes with girth between 10 and 
24. This representation is compact and it is possible to apply 
low-complexity encoding, searching, and decoding procedures 
well developed for convolutional and tailbiting block codes 



[ 32 1, ||33j . In particular, the presented codes with girth 10 
and 12 are shorter than the codes presented in [26| and [8|, 
ll22l . respectively. Moreover, our codes with girth 14 to 18 
are shorter than the corresponding codes presented in [24|. 
Section HXl concludes the paper with some final remarks. 

II. Generator and Parity-Check Matrices 

Consider a rate R — b/c binary convolutional code C with 
the semi-infinite generator matrix 

(Go Gi ... G mg \ 

Go Gi ... G mg (1) 

of memory m g where Gi, i = 0, 1, ... , rn g , are b x c binary 
matrices. Its semi-infinite syndrome former 

H T ■■■ H m 

H T 



H T = 



HI 



HI 



(2) 



of memory m, where in general m ^ ; 
0,1, ... ,m, are (c — b) x c binary matrices, and T denotes 
transpose. Clearly G and H satisfy 



and 



where 



GH 1 



vH 1 







uG 



(3) 



(4) 



(5) 



is the code sequence and u is the information sequence. 

Next we tailbite the semi-infinite generator matrix ([TJ to 
length M c-tuples, where M > maxjrn, m g }. Then we obtain 
the Mb x Mc generator matrix of the quasi-cyclic (QC) block 
code B as 



/ Go 



G T 



G 1 ... 

Go Gi 
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Gm a — 1 G r < 
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(6) 



Gi G2 

Every cyclic shift of a codeword of B by c places modulo 
Mc is a codeword. The corresponding tailbiting parity-check 
matrix is the M(c — 6) x Mc matrix 
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Hi Ho 

: iii 



iijB — 



H„ 



H a 

Hi Ho 



H, n -i 

H r n 



H 2 



H„ 



H„ 



Hi 



H m Hrn — 1 



(7) 
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It is easily shown that Gtb and Hjb satisfy 
Gtb^tb = 



(8) 



given that (O is fulfilled. 

The parity-check matrix for the convolutional code C can 
also be written as the (c — b) x c polynomial matrix 



H(D) = Ha+H^ + ^D 2 
or, equivalently, as 



H{D) 



( h u (D) h 12 (D) 
h 21 {D) h 22 {D) 



H m D" 



hic(D) \ 
h 2c {D) 



(9) 



(10) 



\h {c _ b)1 (D) h [c ^ b)2 (D) ... h [c _ b)c (D) J 



In the sequel we mostly consider parity-check matrices 
with only monomial entries hij(D) = D Wij of degree my, 
where my are nonnegative integers. Clearly, such a parity- 
check matrix H(D) can be represented by its degree matrix 
W = (wij), i = 1, 2, . . . , c — 6 and j = 1, 2, . . . , c. Note that 
starting from Section IVIIII we will relax the restriction to only 
monomial entries and also include zero entries. 



Example 1: Consider the rate R 
with parity-check matrix 



1/4 convolutional code C 



1111 

H(D) = ( 1 1 D D 
1 D 1 D 



(11) 



whose degree matrix is 



W 




11 
10 1 



Tailbiting (fTTT i to length M = 2, we obtain the tailbitten 6 x ; 
parity-check matrix of a QC LDPC block code 
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(12) 



In particular, every cyclic shift of a codeword by c = 4 places 
modulo Mc = 8 is a codeword. 

Due to the restriction to monomial entries in H(D), H TB 
is (J, K)-regular, that is, it has exactly J and K ones in 
each column and row, respectively. Moreover, to fulfill the 
low density criterion, M has to be much larger than J and K, 
and thus the matrix iJxB is sparse. 

Note that the first c columns of ifpB are repeated throughout 
the whole matrix in a cyclicly shifted manner. By reordering 
the columns as 1, c + 1, 2c + 1, . . . , (M - l)c +1, 2, c + 
2, 2c + 2, ... , (M - l)c + 2, etc. and the rows as 1, (c - b) + 
l,2(c-6) + l,...,(M-l)(c-6) + l, 2,(c-6)+2,2(c-6) + 
2, . . . , (M — 1) (c — b) + 2, etc. we obtain a parity-check matrix 



of an equivalent (J, if )-regular LDPC block code constructed 
from circulant matrices 







Iw 12 




H c = 




1-W22 













(13) 



where Wij are the entries of the degree matrix W and I Wi . 
denotes an M x M circulant matrix, that is, an identity matrix 
with its rows shifted cyclically to the left by my positions. 
Note, that the (J, if )-regular LDPC block code determined by 
iic is not quasi-cyclic, although equivalent to the QC block 
code determined by iixB- 

Example Q] (Cont'd): We return to dT2b in Example Q] and 
reorder the columns as 1,5,2,6,3,7,4,8 and the rows as 
1, 4, 2, 5, 3, 6. Then we obtain the equivalent rate R = 1 — 6/8 
(3, 4)-regular LDPC block code with parity-check matrix 
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III. Graphs and Biadjacency Matrices 

A graph Q is determined by a set of vertices V = {vi} and 
a set of edges £ = {e^}, where each edge connects exactly 
two vertices. The degree of a vertex denotes the number of 
edges that are connected to it. If all vertices have the same 
degree c, the degree of the graph is c, or, in other words, the 
graph is c-regular. 

Consider the set of vertices V of a graph partitioned into t 
disjoint subsets Vfe, k = 0, 1, . . . , t— 1. Such a graph is said to 
be t-partite, if no edge connects two vertices from the same 
set Vk, k = 0, l,...,t - 1. 

A path of length L in a graph is an alternating sequence 
of L + 1 vertices Vi, i = 1, 2, ...,£ + 1, and L edges e^, 
i = 1, 2, . . . , L, with ei ^ e%+i- If the first and the final vertex 
coincide, that is, if v\ = u^+i, then we obtain a cycle. A cycle 
is called simple if all its vertices and edges are distinct, except 
for the first and final vertex which coincide. The length of 
the shortest simple cycle is denoted the girth of the graph. 
In |[T8l it was proven that the girth of a graph coincide 
with the minimum distance of the block code, whose parity- 
check matrix corresponds to the incidence matrix of the graph. 
Moreover, the girth determines the number of independent 
iterations in belief-propagation decoding (TJ. 

Every parity-check matrix H of a rate R = k/n LDPC 
block code can be interpreted as the biadjacency matrix [34| 
of a bipartite graph, the so-called Tanner graph [35 1, having 
two disjoint subsets Vo and Vi containing n and n — k vertices, 
respectively. The n vertices in Vo are called symbol nodes, 
while the n — k vertices in V\ are called constraint nodes. If 
the underlying LDPC block code is (J, K )-regular, the symbol 
and constraint nodes have degree J and K, respectively. 
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Si S2 S3 S4 S 5 S6 S7 Sg 




Cl C 2 C 3 C4 C 5 C 6 

Fig. 1. Tanner graph with 8 symbol nodes (Si, i = 1,2, ... ,8) and 6 
constraint nodes (cj, i = 1, 2, . . . , 6). 




Ci C2 C3 



Fig. 2. Bipartite graph with 4 symbol nodes (s^, i = 1,2,3,4) and 3 
constraint nodes (c;, 2 = 1,2,3). Since the edges are labeled according to 
fll6t . this corresponds to a voltage graph. 



Consider the Tanner graph with the biadjacency matrix ifpB, 
corresponding to a QC (J, i"T)-regular LDPC code, obtained 
from the parity-check matrix of a tailbiting LDPC block code. 
Clearly, by letting the tailbiting length M tend to infinity, 
we obtain a convolutional parity-check matrix H(D) ( TTOb of 
the parent convolutional code C. In terms of Tanner graph 
representations, this procedure corresponds to unwrapping the 
underlying graph and extending it in the time domain towards 
infinity. Hereinafter, we will call the girth of this infinite 
Tanner graph the free girth and denote it gf lee . 

Example Q] (Cont'd): Interpreting <TT~4T > as a biadjacency ma- 
trix, we obtain the corresponding Tanner graph Q as illustrated 
in Fig. 03 with 8 symbol nodes and 6 constraint nodes, having 
girth g = 4. 

IV. Base Matrices, Voltages, and their Graphs 

A binary matrix B is called base matrix for a tailbiting 
LDPC block code if its parent convolutional code with parity- 
check matrix H(D) has only monomial or zero entries and 
satisfies 

B = H ( D )\d=i ^ 

that is, all nonzero entries in H(D) are replaced by D° = 1. 
Different tailbiting LDPC block codes can have the same base 
matrix B. 

The base graph Qg follows as the bipartite graph, whose 
biadjacency matrix is given by the base matrix B. Denote 
the girth of such a base graph g%. The terminology "base 
graph" originates from graph theory and is used, for example, 
in ll36l . It differs from the terminology used in [8|, [31], where 
protograph or seed graph are used. 

Consider the additive group (T, +), where T — {7}. From 
the base graph Q B — {£b,Vb} we obtain the voltage graph 
|37|, [38] Qy = {£b,Vb,F} by assigning a voltage value 
j(e ) v, v') to the edge e connecting the vertices v and v', 
satisfying the property j(e, v, v') = — 7(e, v', v). Although 
the graph is not directed, the voltage of the edge depends on 
the direction in which the edge is passed. Finally, define the 
voltage of the path to be the sum of the voltages of its edges. 

Let Q = {£ , V} be a lifted graph obtained from a voltage 
graph Qy, where £ = £s x V and V = Vb x V. Two vertices 
(v,j) and (v',j') are connected in the lifted graph by an 



edge if and only if v and v' are connected in the voltage 
graph Qy with the voltage value of the corresponding edge 
given by 7(e,w,w') = 7 — 7'. It is easy to see that cycles 
in the lifted graph correspond to cycles in the voltage graph 
with zero voltage. Consequently, the girth gv of a voltage 
graph follows as the length of its shortest cycle with voltage 
zero, which is equal to the free girth gf ree O, |23l . A voltage 
assignment corresponds directly to selecting the degrees of the 
parity-check monomials in H(D). 

In the following we start from a base graph Q% and use 
a voltage assignment based on the monomial degrees Wij of 
the degree matrix W to determine the corresponding voltage 
graph Qy. The edge voltage from the constraint node Q to the 
symbol node Sj is denoted by fiij while the opposite direction 
from symbol node Sj to constraint node Ci follows as flji, 
i = 1, 2, . . . , (c — b) and j = 1, 2, . . . , c, where 

A% = wn (i6) 

When searching for LDPC convolutional codes with given 
free girth gf ree , we use integer edge voltages, that is, we 
deal with an infinite additive group. However, when searching 
for QC LDPC block codes with given girth g, obtained by 
tailbiting a parent convolutional code to length M, we use a 
group of modulo M residues, that is, ( fl6] l is replaced by 

Uij — Win mod M 
_ 3 3 (17) 
Hji — —Wij mod M 

The definitions of path and cycle in a voltage graph coincide 
with those in a regular graph, except for the additional 
restriction that two neighboring edges may not connect the 
same nodes in reversed order. The voltage of a path or cycle 
within a voltage graph, follows as the sum of all edge voltages 
involved. 

Example\J] (Cont'd): The bipartite graph whose biadjacency 
matrix is given by the base matrix B of the rate R = 1/4 
(3, 4)-regular LDPC convolutional code C is illustrated in 
Fig. I2 As the edges are labeled according to ([Tol l, Fig. |2] 
corresponds to a voltage graph with girth gy = 4 (for example, 
si — > ci — >• «2 — >• C2 —> Si), The edge from, for example, 
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constraint node c 2 to symbol node S3 is labeled according to 

A*23 = ~M32 = W 2 3 = 1 

As the free girth of the infinite Tanner graph, corresponding 
to the parent convolutional code C, determined by the convo- 
lutional parity-check matrix H{D) is equal to the girth of the 
voltage graph, we can conclude that (?f ree = gy = 4. 

If we neglect all edge labels, we would obtain the corre- 
sponding base graph. 

V. Bounds on the girth and the minimum distance 
of (J > 3, K) QC LDPC block codes 

There are a number of approaches which can be applied 
to construct and search for QC (J = 2, if) -regular LDPC 
block and convolutional codes [12] or the bipartite graphs 
constructed by their incidence matrices. Since every LDPC 
convolutional code can be represented by a bipartite Tanner 
graph using the biadjacency matrix, these techniques can be 
applied to (J > 3, K) QC LDPC codes. Moreover, bounds on 
the girth and the minimum distance of ( J = 2, K) QC LDPC 
codes [12 1 can be generalized to an arbitrary J. 

Theorem 1: The minimum distance G? m j n and the girth g of 
an (n, k, d m \ n ) QC LDPC block code B obtained from a rate 
R = b/c convolutional code C with free distance df ree and 
girth (7f ree by tailbiting to length M are upper-bounded by the 
inequalities 

^min — ^frec 
9 < fffrcc 

Proof: The first statement follows directly from the fact 
that any codeword v(D) of the tailbiting block code B, 
obtained from the parity-check matrix H(D) of the parent 
convolutional code C, satisfies 

v{D)H T {D) = mod {D M - 1) (18) 

Since the parent convolutional code C satisfies ( TT8l without 
reduction modulo (D M — 1) and reducing modulo (D M — 1) 
does not increase the weight of a polynomial, the first state- 
ment follows directly. 

For the second statement we consider the voltage graph Qy 
representation of the parent convolutional code C with girth 
ffv = <?free together with the Tanner graph representation of the 
QC LDPC block code B with girth g. Similar to the relations 
between the free distance df ree and the minimum distance d^, 
there exists a relation between each cycle within the voltage 
graph Qy of the parent convolutional code and the Tanner 
graph Q of the corresponding block code obtained by tailbiting 
to length M, The edge voltages for every cycle in Qy have 
to sum up to zero. Similarly, every cycle in Q corresponds to 
a cycle in Qy such that its edge voltages have to sum up to 
zero modulo M. With the same argument as before it follows 
directly that 

9<9V = fffree 



In lfl2l a lower bound on the girth of a voltage graph gy 
was found via the girth of the corresponding base graph g B for 
ordinary graphs. It is straightforward to generalize this bound: 

Consider a base graph of a QC (J > 3, if) -regular 
LDPC convolutional code with girth g B and let d s denote 
the sth generalized minimum Hamming distance of the linear 
M((J ~ 2)c + b) x JMc block code determined by the parity- 
check matrix which corresponds to the incidence matrix of the 
Tanner graph. In order words, d s corresponds to the number of 
nontrivial (not identically zero) positions of an s-dimensional 
linear subcode. 

Theorem 2: There exist a tailbiting length M and a voltage 
assignment, such that the girth g of the Tanner graph for the 
corresponding TB block code of length N = Mc satisfies the 
inequality 

g > 2 max {g B + \g B /2] ,d 2 }> 3g B (19) 

where d 2 is the second generalized minimum Hamming dis- 
tance, that is, the minimum support of a subcode of dimension 
two. We have equality in ( fT9l ), if the underlying base graph 
consists of two connected cycles, having at least one common 
vertex. 

Proof: According to Theorem Q] any cycle in the Tanner 
graph of a QC LDPC block code corresponds to a cycle of 
the same length in the voltage graph. As the labels of the 
voltage Tanner graph can be freely chosen, it is enough to 
prove that there is no zero cycle shorter than 2(g B + \g-&/2~\), 
that is, no such cycle whose voltage is zero regardless of the 
labeling of the base graph. In particular, such a cycle is also 
known as an inevitable cycle Q or balanced cycle [8|. The 
number of times each edge in such a cycle of the voltage 
graph is passed in different directions has to be even. This 
cycle cannot be simple, since in a simple cycle each edge is 
passed in one direction only. Hence, the cycle passes through 
the vertices of a subgraph which contains at least two different 
cycles, corresponding to two different nonzero codewords. The 
minimum distance of the base graph is g B , and, according to 
the Griesmer bound, the smallest length of a linear code with 
two nonzero codewords of minimum distance d is d+ \d/2~], 
and, hence, the first lower bound of inequality ( fT9l ) follows. 

Consider the second lower bound. The definition of the 
second generalized minimum Hamming distance implies that 
the smallest subgraph with two cycles has to have at least d 2 
edges. Thus, the second of the two lower bounds gives the 
precise value of the girth of a subgraph containing two con- 
nected cycles, having at least one common vertex. Otherwise, 
d 2 is a lower bound. ■ 

The bounds are tighter than the 3i?b bound [7], [36| but not 
tight if the shortest non-simple cycle consists of two simple 
cycles connected by a path. 

Finally, we want to recall an upper bound on the achievable 
girth and minimum distance. We start by reformulating the 
theorem on the achievable girth by Fossorier lfl8l and thereby 
generalize it to include base matrices with zero elements. 

Theorem 3: Consider the parity-check matrix H(D) of a rate 
R = b/c convolutional code with base matrix B. Denote the 
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corresponding base graph Qg and let B' be the 2x3 submatrix 

B '={ \ I I) < 20 > 

If the base matrix B, after possibly reordering its rows and 
columns, contains the submatrix B', then the girth gy of the 
corresponding voltage graph Qy is upper-bounded by 

g Y < 12 (21) 

regardless of the voltage assignment. 

Proof: The subgraph determined by the 2x3 submatrix 
B' contains 3 symbol nodes, 2 constraint nodes, and 6 edges. 
Moreover, there exist 3 shortest cycles of length 4. Thus, the 
base graph Qb has girth <?b = 4 and its second generalized 
Hamming distance is di = 6. Applying Theorem[2] we obtain 
the precise value of the achievable girth as 2d-2 = 12, which 
completes the proof. ■ 

For parity-check matrices with only nonzero monomial 
entries, the inequality ( TSTT i was proven in |fl~8). 

Moreover, let H(D) be the parity-check matrix of a rate 
R = b/c (J, if) -regular LDPC convolutional code with free 
distance df Ke . By tailbiting to length M we obtain a QC LDPC 
block code of block length M c and minimum distance d m i n . As 
proven in |2| for parity-check matrices without zero elements 
and reformulated in [11 1 for parity-check matrices with zero 
elements, the corresponding minimum distance d m \ B can be 
upper-bounded by 

drain < dfree < (c-6+1)! (22) 

For parity-check matrices with only nonzero monomials, the 
inequality simplifies to (J + 1)!. 

VI. Searching for QC LDPC block codes with 

LARGE GIRTH 

When searching for QC LDPC block codes with large girth, 
we start from a base graph of a rate R — b/c (J, if) -regular 
LDPC convolutional code. Using the following algorithm, we 
determine a suitable voltage assignment based on the group 
of nonnegative integers, such that the girth of this voltage 
graph is greater than or equal to a given girth g. Afterwards 
we replace all edge labels by their corresponding modulo M 
residuals, where we try to minimize M while preserving the 
girth g. Using the concept of biadjacency matrices leads to 
the corresponding degree matrix W and we obtain the parity- 
check matrix of a convolutional code whose bipartite graph has 
girth g = gf lee . Tailbiting to lengths M, yields the rate R = 
Mb/Mc QC LDPC block code whose parity-check matrix is 
equal to the biadjacency matrix of a bipartite graph with girth 
9- 

The algorithm for determining a suitable voltage assignment 
for a base graph consists of the following two main steps: 

1) Construct a list containing all inequalities describing 
cycles of length smaller than g within the base graph. 

2) Search for such a voltage assignment of the base graph 
that all inequalities are satisfied. 

The efficiency of the second step, searching for a suitable 
voltage assignment, depends on the chosen representation for 
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Fig. 3. A bipartite voltage graph with 4 symbol nodes (sj, i = 1,2,3,8) 
and 3 constraint nodes (c;, i = 1, 2, 3) with its edges labeled according to 

the list of inequalities determined during the first step. In 
general, when searching for all cycles of length g roughly 
(J — l) 9 different paths have been considered. However, by 
using a similar idea as in l33l when searching for a path within 
a trellis, we create a tree of maximum depth g/2 and search 
only for identical nodes within the tree and thereby reduce the 
complexity to roughly (J — I) 9 / 2 . 

Creating a tree structure 

Utilizing the base graph of a rate R — b/c (J, if) -regular 
LDPC convolutional code, with c symbol and c — b constraint 
nodes, we construct a separate subtree starting with each of 
the c symbol nodes. 

Before describing the algorithm, we have to introduce some 
notations. A node in the tree will be denoted by £ and has a 
unique parent node £ p . The underlying base graph is bipartite, 
that is, every node £ in the tree with £ G V, is only connected 
to nodes £' G Vj with i,j G {0, 1}, i ^ j, where Vo and Vi 
are the sets of symbol and constraint nodes, respectively. In 
other words, a symbol node is only connected to constraint 
nodes and vice versa. Moreover, every node £ is characterized 
by its depth £(£) and its number ra(£), where n(£) = i follows 
directly from £ = Sj or £ = a depending on whether its depth 
£(£) is even or odd. In particular, £ G Vt(£) mod 2- 

Having introduced those basic notations, we can grow c 
separate subtrees, with the root node £i yl - oot of the ith subtree 
being initialized with £ G Vo and depth £(£) = 0. Extend 
every node £ G V at depth £(£) = n with i = n mod 2 by 
connecting it with the nodes £' G Vj+i mod 2 at depth n + 1 
according to the underlying base graph, except £ p which is 
already connected to £ at depth n — 1. 

Finally we label the edges according to ( TToT l and obtain the 
voltage for node £ in the ith subtree as the sum of the edge 
voltages of the path £;. r oot — > £,■ 

Clearly, all subtrees together contain all paths of a given 
length in the voltage graph. Moreover, taking into account 
that the girth g of bipartite graphs is always even, we can 
conclude that in order to check all possible cycles of length 
at most g — 2 in the voltage graph, it is sufficient to grow 
the corresponding c subtrees up to depth (g — 2)/2 and to 
construct voltage inequalities for all node pairs (£,£') in the 
same subtree i with the same number n(£) = n(£') and depth 
£(£) = £(£') but different parent nodes ^ ^ £' p - 

The corresponding voltage inequality for the node pair 
(£,£') follows directly as the difference between the voltages 
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Fig. 4. A tree representation with maximum depth two, starting with symbol 
node si. 



for the paths from ^., oot to £ and respectively, that is, 

(£i,root £) — (£i,root — ► £,')■ 

If there exists a cycle g' < g, then at depth g'/2 there exists 
at least one pair of nodes (£, £'), whose corresponding voltage 
inequality is not satisfied, that is, is equal to zero. If there is 
no cycle shorter than g in the voltage graph, then there are no 
such pairs, and all voltage inequalities are satisfied. 

Example 2: Consider the rate R — 1/4 (3, 4)-regular LDPC 
convolutional code given by (fTTb . The corresponding base 
graph, with four symbol nodes Si £ Vo, % = 1, 2, 3, 4, and three 
constraint nodes Cj 6 Vi, i = 1,2, 3, is illustrated in Fig. [3] In 
the following, we shall search for a set of edge labels, that is, 
monomial degrees in W, such that the corresponding voltage 
graph has at least a given girth g. Thus, we label the branches 
by the general edge voltages according to (TToT i and obtain a 
bipartite voltage graph. 

In order to find a suitable labeling for the edge voltages 
from the ith constraint node Cj to the jth symbol node Sj, that 
is Hij, i = 1,2,3, j = 1,2,3,4, whose underlying voltage 
graph has at least girth g = 6, we have to grow 4 subtrees up 
to length (g — 2)/2 = 2, with their root nodes being initialized 
by Si, i = 1,2,3,4. 

The subtree with root node si is illustrated in Fig. |4] Clearly, 
at depth £(£) = 1 there are no identical nodes, while at depth 
£(£) = 2 there are 3 x =9 pairs of identical nodes (n(£) = 
but with different parents. Taking into account that a 
similar subtree is constructed using the remaining three symbol 
nodes S2, S3 and S4 as root nodes, we obtain in total 36 = 4 x 9 
node pairs, which all correspond to a voltage inequality. 

For example, the voltage inequality obtained by checking 
all node pairs (£,£') with £ = S2, that is, n(£) = 2, at depth 



£(£) = 2 in the subtree starting with symbol node si, are 

-Mil + M12 - M22 + M21 
-fin + M12 - A*32 + M31 7^ 
-M21 + M22 - M32 + A*3i 

Note that amongst all 36 voltage inequalities, there are only 
18 unique ones. 

Algorithm TR: Constructing a tree representation. 

1) Grow c separate subtrees according to the underlying 
base graph up to depth g/2 — 1, with the root node 
£i,root °f me *th subtree being initialized with £ G Vo 
and depth £(£) = 0. 

2) Extend every node £ e V,; at depth £(£) = n with i = n 
mod 2 by connecting it with the nodes £' € Vj+i m od 2 
at depth n + 1 according to the underlying base graph, 
except £ p which is already connected to £ at depth n—1. 
Denote the set of all nodes within the ith subtree by %. 

Searching for a suitable voltage assigment 

Using the obtained subtrees 71, i = 1,2, ... ,c, with maxi- 
mum depth g/2 — 1, we have found all cycles of length smaller 
than or equal to g — 2 as well as their corresponding voltage 
inequalities. 

The same cycle might be found several times within the c 
subtrees. Moreover, two different cycles can correspond to the 
same voltage inequality. 

We continue by creating a reduced list C of node pairs 
(£>£') °f a U c subtrees %, i — l,2,...,c, containing all 
unique voltage inequalities. Thereby, we remove all duplicate 
cycles, as well as different cycles corresponding to the same 
voltage inequality. Using the reduced list L we can reduce 
the obtained c subtrees %, i = 1, 2, . . . , c, in a similar way by 
removing all nodes, not participating in any of the cycles in C, 
and denote the reduced subtree by T^mm- In other words, we 
remove all nodes in 77, i = 1, 2, . . . , c, which only participate 
in already known cycles or new cycles with already known 
voltage inequalities. 

In the following we present two different approaches for 
finding suitable edge labels (edge voltages), which we shall 
denote as Algorithm A and Algorithm B. 

In Algorithm A, we label the edges of the reduced subtrees 
7i,min, i = 1,2, . . . , c, with a set of randomly chosen voltages. 
For every node pair (£,£') in the list C, we calculate the 
voltage of the corresponding cycle as the difference of the 
path voltages Ci.mot — > £ and £;. root — > If none of these 
cycle voltages is equal to zero, the girth of the underlying 
base graph with such a voltage assignment is greater than or 
equal to g. 

In Algorithm B, we discard the list L and only focus on 
the c reduced subtrees 7i, m in- After labeling their edges with a 
set of randomly chosen voltages, we sort the nodes £ of each 
subtree according to their path voltage £i !roo t — > £■ If there 
exists no pair of nodes (£, £') with the same path voltage, 
number n(£) = n(£')> and depth £(£) = £(£'), but different 
parent nodes £ p ^ £' p , the girth of the underlying base graph 
with such a voltage assignment is greater than or equal to g. 
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A formal description of those two algorithms is given below: 

Algorithm A: Constructing a system of voltage inequalities and 
searching for an optimum voltage assignment using a list. 

1) Create a reduced list C of node pairs (£,£') for all c 
subtrees %, i — 1,2, ...,c, containing all node pairs 
(£,£') with a unique voltage inequality, having the same 
number n(£) = n(f'), depth £(£) = but different 
parent nodes £ p ^ £' p . 

2) Reduce each of the c subtrees % by removing all 
nodes, which do not participate in any of the found 
cycles corresponding to the voltage inequalities in C, 
and denote the reduced subtree structure by Ti^m- 

3) Assign randomly chosen voltages to the edges of all 
trees and perform the following steps: 

a) Find the voltages for all paths leading from the root 
node £i, root of the zth reduced subtree %,mm to all 
nodes £ 6 7i,min, i = 1,2,.. .,c. 

b) Determine the voltage inequality for all cycles 
(£, f) G C, as the difference of the corresponding 
path voltages in 7i jm i n , i — 1, 2, . . . , c, computed 
previously. 

c) If all voltage inequalities are satisfied, the girth 
of the underlying base graph with such a voltage 
assignment is greater than or equal to g. 

Algorithm B: Constructing a system of voltage inequalities and 
searching for an optimum voltage assignment using a tree. 

1) Construct the reduced list C and the reduced subtrees 
7i,min, i = 1, 2, . . . , c, as in Algorithm A without storing 
the corresponding list C. 

2) Assign randomly chosen voltages to the edges of all 
trees and perform the following steps: 

a) Find the voltages for all paths from the root node 
£j.root to all nodes within 7I. m in, i = 1, 2, . . . , c, and 
sort all elements within 77,min according to their 
voltages. 

b) Search for a pair of nodes (£, £') in the sorted list 
with the same path voltage, number n(£) = 

and depth = £(£,'), but different parent nodes 

c) If no such pair exists, then the girth of the corre- 
sponding voltage graph with such a voltage assign- 
ment is greater than or equal to g. 

Complexity 

Denote the sum of all nodes in the reduced tree 7i, m in, i = 
1,2, ... ,c, and the number of unique inequalities in the list C 
by Nj and JVl, respectively, that is, 

c 

Nj = Y,\n min \ and N L = \£\ 

i=l 

where \X\ denotes the number of entries in the set X. 

Algorithm A requires Nj summations for computing the 
path voltages and JVl comparisons for finding cycles, leading 
to the complexity estimate Nj + Nl- Algorithm B requires 
the same number of JVt summations for computing the path 



TABLE I 

Complexity of a search for suitable voltage assignment for 
qc ldpc block codes with girth g < 12 



K 


a 


= 8 


- 


= 10 


a 


= 12 


Nj 




JVt 


N L 


Nj 


N L 


4 


53 


42 


150 


231 


269 


519 


5 


93 


90 


286 


645 


581 


1905 


6 


142 


165 


485 


1470 


1060 


5430 


7 


200 


273 


759 


2919 


1742 


12999 


8 


267 


420 


1120 


5250 


2663 


27426 


9 


343 


612 


1580 


8766 


3859 


52614 


10 


428 


855 


2151 


13815 


5358 


93735 


11 


522 


1155 


2845 


20790 


7210 


157410 


12 


625 


1518 


3674 


30129 


9446 


251889 



voltages, roughly Nj log 2 Nj operations for sorting the set, 
and JVj comparisons, leading to a total complexity estimate 

of JVt log 2 JVt. 

In Table|I]the values of JVt and JVl are given when searching 
for suitable voltage assignment for a (J, A")-regular rate R = 
1 — J/K QC LDPC convolutional codes with J = 3 and 
arbitrary K > 4 and girth g constructed from all-ones base 
matrices. In this case, up to g = 10, Algorithm A is preferable, 
while when searching for voltage assignment with girth g > 
12, Algorithm B should be used. 

In the general case we have to consider all node pairs, and 
as JVl is roughly JV| we conclude that Algorithm B performs 
asymptotically better (when JVt — > oo). 

VII. Minimum Distance of QC LDPC Codes 

Usually the girth of the Tanner graph of an QC LDPC block 
code is considered to be the most important parameter that 
affects the performance of belief-propagation decoding, as it 
determines the number of independent iterations (TJ. There- 
fore, most research is focused on finding QC LDPC block 
codes with large girth, while their corresponding minimum 
distance is mostly unknown. In [14] it was shown, that the 
performance of belief -propagation decoding algorithms at high 
SNRs depends on the structure and the size of the smallest 
absorbing sets, which however can be upper-bounded by the 
minimum distance. This is the rationale for computing the 
minimum distance of the shortest known QC LDPC block 
codes as well as our search for QC LDPC codes with both 
large girth and large minimum distance. 

Our method of calculating the minimum distance is based 
on the well-known fact that the minimum distance of a linear 
block code B with parity-check matrix H is equal to the 
minimum number of columns of H which sum up to zero. 

Consider the M(c — b) x Mc parity-check matrix 7?tb of 
the (J, i^Q-regular rate R = Mb/Mc tailbiting block code 
B with block-length JV = Mc ©■ Starting with each of the 
first c columns of ifpB as a root, we will construct c separate 
trees, where each node £ is characterized by its depth and 
partial syndrome state column-vector <r(£). 
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TABLE II 

Degree matrices for QC LDPC codes with girth g = 6 



K | (n, fc) | d min | M | W_ 

short codes 



4 


(20, 7) 


6 


5 


12 4 
3,1,2 


5 


(25, 12) 


6 


5 


1.234 
3,1,4,2 


6 


(42, 23) 


4 


7 


1 2 3 4 6 
3,5,2,1,4 


7 


(49, 30) 


4 


7 


1 2 3 4 5 6 
3,5,2,1,6,4 


8 


(72,47) 


4 


9 


1 2 3 4 5 7 8 
3,6,2,1,8,5,4 


9 


(81,56) 


4 


9 


1,2,3,4,5,6,7,8 
3,6,2,1,8,7,5,4 


10 


(110,79) 


6 


11 


1,2,3,4,5,6,8,9, 10 
3,1,7,2,10,9,4,6,5 


11 


(121,90) 


4 


11 


1,2,3,4,5,6,7,8,9,10 
3,1,7,2,10,9,8,4,6,5 


12 


(156, 119) 


6 


13 


1,2,3,4,5,6,7,8, 10,11,12 
3,1,8,2,9, 12,4,11,5,7,6, 


large distance codes 


4 


(92, 25) 


22 


23 


1,2,4 
5,3, 12 


5 


(245, 100) 


22 


49 


1,3, 10, 14 
40,31,33,30 


6 


(414, 209) 


22 


69 


3,4,21,26,67 
34, 15, 64, 33, 44 


7 


(763, 438) 


22 


109 


1,3,11,15,45,93 
101,34, 18,9,1,4 


8 


(1224, 767) 


22 


153 


2,10,26,57,89,4,49 
22,19,5,23,61,90, 123 



Initialize the partial syndrome state of the root £i jr0 ot of 
the ith tree with column i of the corresponding parity-check 
matrix, that is, <r(^i, loot ) = hi, i = 1,2, ...,c. Then grow 
each tree in such a way, that every branch between any two 
nodes £ and £' is labeled by a column vector hj, j ^ i, such 
that <t(£') = + hj, where every branch label on the path 
£i,root C does not occur more than once. 

Consider now a certain node £ with nonzero state <x(£) = 
(o-i(0 0-2(0 ■ • ■ <T( C - 6 )(£)) T , where Cj (£), j = 1,2, . . . , c-6 
is aMxl column vector. If we assume that the fcth position 
of the column-vector <r(£) is nonzero, then there are at most 
K — 1 columns which can cancel this nonzero position and 
have not been considered previously. Therefore, every node £ 
has at most K — 1 children nodes per nonzero position. 

Such a tree would grow until all possible linear com- 
binations have been found. Therefore, we assume that the 
minimum distance is restricted by d m [ n < t, that is, the 
maximum depth of the tree is t — 1. Consequently, a node 
£ at depth will not be extended, if the number of nonzero 
positions in its partial syndrome state column-vector <x(£) 
exceeds J(t — £(£) — 1), since at most J ones can be canceled 
by each branch. 

By initially reordering the columns of the parity-check 
matrix i/jB in such a way that each of the c— b nonoverlapping 
blocks of M rows contains not more than a single one per 



column, we can strengthen the stopping criterion as follows: 
A node £ at depth £(£) will not be extended, if the number of 
nonzero positions in each of its partial syndrome state column- 
vectors (Tj(^), j = 1,2, ... ,c — b exceeds (t — £(£) — 1), since 
at most 1 one in each block can be canceled by each branch. 
In particular, such a reordering of the parity-check matrix ifpB 
corresponds directly to the parity-check matrix He (TT~3T > of the 
equivalent ( J, K) -regular LDPC block code constructed from 
circulant matrices. 

Algorithm MD: Determine the minimum distance of a rate 
R = b/c (J, K)-regular LDPC block code. 

1) Assume a suitable restriction t on the minimum distance 

^min ^ t' 

2) Grow c separate trees as follows: 

a) Initialize the root node of the ith tree by 
f(£ r0 ot,i) = hi with depth £(£) = 0. 

b) Extend all nodes £ as long as the Hamming weights 
of their partial syndrome states Wn(<r(£)) < J(t — 
£(£) — 1) (Note, for codes with blocks of M rows 
containing only a single one, this criterion can be 
strengthen to w H (er .,(£)) < (t - £(£) - 1), j = 
l,2,...,c-b). 

c) The minimum distance e£ m j n follows directly as 

d min = nun {^(0 | <t(0 = 0} 

If there is no node £ whose partial syndrome state 
<t(£) = 0, then the minimum distance is lower- 
bounded by d m i„ > t. 

VIII. Search results 

When presenting our search results for QC ( J = 3, ir- 
regular LDPC block codes with different girth we will distin- 
guish two cases. 

We started by searching for QC ( J = 3, if)-regular LDPC 
block codes using an all-one base matrix B, applied the 
algorithms as described above and obtained QC (3, iQ-regular 
LDPC block codes with girth g = 6, 8, 10, and 12 as given in 
Tables iHl-IVl These codes correspond to a parity-check matrix 
H(D) of a convolutional code with only monomial entries as 
given, for example, in (fTOb . 

However, according to Theorem [3] the achievable girth g of 
a QC (J, i"Q-regular LDPC code, constructed in such a way, 
is limited by g < 12. Thus, in order to find QC (J = 3, ir- 
regular LDPC block codes with girth g > 12 as presented 
in Tables [VI] and IVHI we have to allow zero entries in our 
base matrix B. This is a straight-forward generalization of the 
restriction to only monomial entries in the parity-check matrix 
H(D) of the underlying convolutional code. 

Case I: monomial entries 

In Tables HI1- IVl parity-check matrices of short known QC 
( J = 3, irregular LDPC block codes with girth 5 = 6,8, 10, 
and 12 together with those of large minimum distance are 
presented. When searching for such codes, we applied the 
following restrictions to reduce the number of possible voltage 
assignments: 



TABLE III 

Degree matrices for QC LDPC codes with girth g = 8 



K | (MQ | o! mm | M | W_ 

short codes 



4 


(36, 11) 


6 


9 


1,4, 6 
5,2,3 


5 


(65, 28) 
((75,32) [24]) 


10 


13 


1.3, 7, 11 
10,4,5,6 


6 


(108, 56) 
((156, 80) E3) 


10 


18 


2, 3, 5, 7, 9 
4,6, 13, 1, 16 


7 


(147, 86) 


10 


21 


2, 3, 8, 15, 17, 20 
4,6,7,9, 12,13 


8 


(200, 127) 


8 


25 


1,3, 4, 10, 14, 15, 19 
5,6,11,24,2,9, 12 


9 


(270, 182) 


8 


30 


1,3,10, 16,23,25,26,28 
2,6,5,9,8,12, 14,22 


10 


(350, 247) 


8 


35 


2,6,7, 18,19,26,29,31,34 
4,5,3,13,10, 16,12,11,23 


11 


(451,330) 


8 


41 


1, 4, 8, 20, 27, 28, 29, 33, 39, 40 
5,7,6,9, 10, 19,13,21,14,35 


12 


(564, 425) 


8 


47 


3, 7, 8, 22, 24, 27, 29, 35, 40, 41, 43 
6,2,4,5, 14, 16,1,21,28,9,34 


large distance codes 


4 


(116, 31) 


24 


29 


3,14,21 
7,1,17 


5 


(225,92) 


24 


45 


1,3, 10, 14 
40,31,33,30 


6 


(431,218) 


24 


72 


3,4,21,26,67 
34, 15,64,33,44 


7 


(777, 446) 


24 


111 


3,11,15,45,93,110 
34, 18,9, 1,4, 101 


8 


(1280,802) 


24 


160 


2,4, 10,26,49,57, 89 
22,90, 19,5, 123,23,61 


9 


(1386,926) 


20 


154 


6,9,26,65,79,99, 124, 153 
24, 16,14,1,46,62,137, 84 



• As the girth of a voltage graph is defined as the shortest 
cycle with voltage zero, and the sign of the voltage 
depends on the direction of the edge, we can add the same 
arbitrary offset to the voltage of all edges being connected 
to the same node. Thus, without loss of generality, we 
set the voltage of all edges connected to one specific 
symbol node as well as all edges connected to one specific 
constraint node to voltage zero. (For consistency with 
codes constructed from Steiner Triple System, that will 
be introduced later, we choose the first symbol node and 
the last constraint node. This corresponds directly to a 
degree matrix with zeros in its first column and last row.) 
For example, the degree matrix of the (J = 3,K = 4) 
QC LDPC block code with girth 5 = 8 from Table EI] 
follows directly as 

/ 1 4 6 \ 
W= 5 2 3 
\ J 

• Furthermore, to reduce the number of only permuted 
degree matrices, we assume that 

- The first row is sorted in ascending order. 

- When sorting the first and the second row in ascend- 
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TABLE IV 

Degree matrices for QC LDPC codes with girth g = 10 



K I (M) |°U,| M I W_ 

short codes 



4 


(148, 39) 
((144,38) El) 


14 


37 
(39 (26)) 


I, 14, 17 

II, 6,2 


5 


(305, 124) 
((550,222) [24]) 


24 


61 

(6i ED) 


2, 20, 54, 60 
26,16,31,48 


6 


(606, 305) 
((780,392) (24l) 


24 


101 

(103 HO) 


2,24,25,54,85 
21,15,11,8,59 


7 


(1113,638) 


24 


159 
(160 (2D) 


2 14 27 67 97 130 
21,24,1,6,75,58 


8 


(1752, 1097) 


24 


219 

(233 HO) 


3, 14,26,63,96,128,183 
24,6,19,46,4,77, 107 


9 


(2871, 1916) 


24 


319 
(329 m) 


ft Q 2fi fi5 QQ 1 5^ 9^ 978 

24,16, 14,1,62,84,200, 137 


10 


(4300, 2912) 


24 


430 


9,11,26,67, 101,161,... 
233,302,395 

9'^ 5 1 54 33 Qfi 
iO, u, J-, o*±x oo, yu, . . . 

120, 104, 244 


11 


(6160, 4482) 


24 


560 
(577 HO) 


2,11,25,62,101,162,225,... 

268,421,492 
24,21,5,55,6,59,178,... 

132,204,311 


12 


(8844, 6635) 




737 
(758 [26]) 


2,22,23,63,101,147,219,... 

322,412,569,601 
16,9,6,58,34,91,126,... 

155,185,298,232 


large distance codes 


4 


(176, 46) 


24 


44 


I, 14, 17 

II, 6,2 



ing order, the second row is lexicographically less 
than the first row. 
- The maximum degree is less than the tailbiting length 
M for which there exists a QC (J = 3, K )-regular 
LDPC block code with given girth g. 
■ QC (J = 3,K = 4)-regular LDPC block codes were 
found by exhaustive search over the previously defined 
set of restricted edge voltages. 
. QC (J = 3, K = AO-regular LDp C block codes with 
A > 4 were obtained by adding one additionally, 
randomly chosen column to the best degree matrices of 
codes with K = A — 1 having the same girth g. The 
maximum degree in this additional column is limited by 
twice the maximum degree of the previous code. 
Using these restrictions, the obtained QC (J = 3, AT)-regular 
LDPC block codes with girth g = 6,8,10, and 12 are 
presented in Tables iHl-IVl 

The first column K denotes the number of nonzero elements 
per row, which corresponds to the number of columns in H(D) 
and W, due to the all-ones base matrix B. As all entries in 
the first column and the last row of the degree matrix W are 
zero, they are omitted in the submatrix W' which is given in 
the corresponding last column. 

Consider now the parity-check matrix H(D) of the rate R = 
1 — J/K convolutional code C, with only monomial entries 
corresponding to the degree matrix W. By tailbiting the semi- 
infinite parity-check matrix H to length M (given in the forth 
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TABLE V 

Degree matrices for QC LDPC codes with girth g = 12 



K\ (ryfc) | M | W_ 

short codes 



4 


(292, 75) 
((444, 113) [24]) 


24 


73 
(97 @) 


2, 25, 33 
18,6,5 


5 


(815, 328) 
((1700,682) E2) 


24 


163 

(i8i ED) 


5, 33, 42, 117 
36,35,25,57 


6 


(1860, 932) 
((4680, 2342) l24l) 


24 


310 
(393 l22l) 


1,24,38, 145, 246 
16,36,5,82,110 


6 


(1836,920) 


24 


306 
(393 [22]) 


9 36 38 1 54 204 
33, 1, 13,54, 123 


7 


(3962, 2266) 




566 
(881 E) 


3,10,33,147, 297, 442 
31,22,4,93,133,219 


8 


(6784, 4242) 




848 
(1493 0) 


A 9zL ^1 '\A 1 \ ^0^ ZLQS fi^9 

^^t, 0±, -L'iO, OUO, 'ic'O, Udi 

32,9,6,70, 130, 193,222 


9 


(12384, 8258) 




1376 

1 zuo I |o| 1 


4,20,32, 160, 284, . . . 
569,794, 1133 

■30 7 1 Q9 1 fiQ 

350,437,645 


10 


(21030, 14723) 




2103 


6,13,28,150, 291,565,... 

678, 1258, 1600 
30, 16, 5,64, 225, 207, . . . 

491,838, 746 


11 


(34507, 25098) 




3137 


9,11,24, 150,306,508,... 

666, 1279, 1765, 1958 
31,28,1,83,131,160,... 

429,550,956,1391 


12 


(56760, 42572) 




4730 


3, 15, 22, 140, 286, 537, . . . 

811,1113, 1878,2524,3349 
31,26,1,66,95,210,373,... 

729,878, 1365,1644 



column), we obtain the parity-check matrix ifpB of an (n, k) 
block code B with minimum distance d m ; n , where (n, k) and 
dmin follow from the second and third column, respectively. 
Note that due to linear dependent rows in i/xB the rank of B 
might be less than M(c — b). 

The codes presented in Tables [II] and [Til] coincide with the 
QC LDPC block codes found by the "hill-climbing" algorithm 
[26 1, but we determined their minimum distance with our 
algorithm described in Section I VIII Tables |IV] and [V] contain 
new QC (J = 3, K )-regular LDPC block codes, which, to 
the best of our knowledge, are shorter than the previously 
known codes obtained from an all-ones base matrix [8|, |19|, 
[22 1, [26 1. In particular, these codes are significantly shorter 
than those presented in ll24l . which are obtained from base 
matrices with zeros. However, due to the zeros in their base 
matrix, the minimum distance of the LDPC block codes in 
ll24l can exceed (J + 1)!. For example, we determined the 
minimum distance of the (444, 113) QC (3, 4)-regular LDPC 
block code with girth g = 12 in [24 1 to be d m i n = 28, while 
the corresponding code in Table [Vl that is, the (292,75) QC 
(3, 4)-regular LDPC block code, has only minimum distance 
d min = 24, but shorter block length. Using the BEAST 11331 . 
we calculated the free distance of the corresponding parent 
convolutional code for the code in [24 \ to be df ree = 46. 
Therefore, using our approach and a larger tailbiting length 
it would be possible to construct corresponding QC (3, 4)- 



regular LDPC block codes with minimum distance up to 46. 

Case 2: monomial or zero entries 

In order to find QC ( J = 3, K ) -regular LDPC block codes 
with girth g > 14, we have to allow zero entries in our base 
matrix B; that is, relax the restriction from only monomial 
entries in H(D) to include also zero entries. According to 
Theorem |2] a code with girth g exists if the corresponding 
base graph has girth g% satisfying (fl9] >. Additionally, as we 
are searching for codes with short block length, we consider 
the shortest possible base matrices B. 

Case 2-1: Steiner Triple Systems 

When searching for QC (J = 3, K ) -regular LDPC block 
codes with girth g — 14, 16, and 18, we started with a 
(shortened) base graph constructed by using Steiner triple 
systems of order n, that is, STS(n) Hg), El, ED. 

For all n, where n mod 6 is equal to 1 or 3, there exists 
a Steiner triple system of order n. Then we construct a 
(J, if) -regular, (c — b) x c base matrix B with entries bij, 
i = 1, 2, . . . , c — b and j = 1, 2, . . . , c, in such a way that the 
positions of the nonzero entries in each column correspond to 
a Steiner triple system of order (c—b). Denote such a (c—b)xc 
base matrix 5 S ts( c -&)- 

Using the obtained (J, K )-regular (c — b) x c base matrix B, 
we search for a set of edge labels, such that the corresponding 
voltage graph has at least girth g. 

In general, it is possible, without loss of generality, to label 
a certain subset of edges of the voltage graph simultaneously 
with zero voltage and thereby decreasing the number of possi- 
ble labelings. The following algorithm constructs a (c — b)xc 
base matrix B based on STS(c — b) and reorders the matrix to 
maximize the number of zero entries in its lower left corner. 
Using such a base matrix, it is always possible to label the last 
nonzero entry in each column with degree zero. Moreover, in 
each of the remaining rows at the top of the base matrix, 
we can label at least one nonzero entry with degree zero. 
(Hereinafter we will always choose at least the first element 
in the remaining rows to be labeled with zero voltage). 

Algorithm STS: Construction of a (J, K)-regular (c — b) x c 
base graph B obtained from STS(c — b). 

1) Initialize a counter u to zero and denote the current row 
and column by s and t, respectively, starting from the 
right-most entry in the last row, that is, s = c — b and 
t = c. 

2) Set the K — u elements in row s and column t, t — 
1, . . . , i — K + u + 1 to one, that is, bij = 1 for i = s 
and j = t,t - 1, . . . , t - K + u + 1. 

3) Choose the remaining J — 1 nonzero positions in each 
of those K — u columns to fulfill the properties of a 
Steiner Triple System. If possible, choose the positions 
b^ to minimize i. In other words, try to avoid using the 
lowest rows s — 1, s — 2, . . ., if possible, despite of the 
restrictions imposed by the Steiner Triple System. 

4) Finally, decrease t by K — u, set s to s — 1, denote the 
number of nonzero elements in the new row s by K — u 



TABLE VI 

Degree matrices for QC LDPC codes with girth g = 14 to 18 



K 


9 


(n,fc) 


M 


Base graph 


w 


1 


14 


(1812,453) 
((2208, 732) 03) 


151 
(184 (24)) 


STS(9) 
(9 X 12) 


0,123, 36,3,2,79,4,7,52,4,1 
0,96,23,11,1,37,12,2,61,1,4 


5 


14 


(9720, 3888) 
((11525,4612) (24]) 


486 


S-STS(13) 
(12 X 20) 


423, 0, 437, 5, 237, 235, 170, 333, 260, 109, 241, 2, 114, 5, 2, 428, 92, 228, 299 
0, 0, 0, 445, 465, 51, 440, 22, 111, 307, 433, 4, 285, 2,1,4, 113, 282, 5 


6 


14 


(29978, 14989) 
((37154, 18579) El) 


1153 
(1429 f24l) 


STS(13) 
(13 x 26) 


1037, 0, 1051, 1105, 933, 1027, 962, 1000, 665, 805, 646, 2, . . . 

906, 5, 2, 1095, 788, 633, 913, 264, 51, 772, 672, 686, 737 
0, 0, 0, 1112, 1132, 51, 1107, 22, 807, 921, 1100, 4, 952, 2, . . . 

1, 4, 905, 949, 5, 0, 1111, 922, 620, 351, 140 


7 


14 


n = 80000000 


800000 


STS(25) 


available at I139\I 


4 


16 


(7980, 1995) 


665 


STS(9) 
(9 x 12) 


0, 468, 99, 3, 2, 305, 43, 9, 251, 3, 2 
0,351,41,6,8,215,18,1,79,1,8 


5 


16 


(51240, 20496) 
((62500,25002) El) 


2562 


S-STS(13) 
(12 x 20) 


937, 0, 1551, 1264, 1670, 2119, 1973, 1960, 1848, 1223, 1806, . . . 

15, 1761, 1, 2, 2175, 1169, 1768, 548 
0, 0, 0, 2367, 2491, 126, 2296, 66, 1197, 582, 2200, 9, . . . 

1836,2,1,0,1757, 1833,4 


6 


16 


(227032, 113516) 
((229476, 114740) [24]) 


8732 
(8826 l24l) 


STS(13) 
(13 x 26) 


8328, 0, 8393, 8106, 7840, 8289, 8143, 8130, 6821, 7393, 6779, 15, 7931, . . . 

1, 2, 8345, 7339, 6741, 7390, 1557, 498, 6357, 5666, 5001, 1684 
0, 0, 0, 8537, 8661, 126, 8466, 66, 7367, 7424, 8370, 9, 8006, 2,1,0,... 

7927, 8003, 4, 0, 8412, 5799, 4553, 2142, 6293 


4 


18 


(32676, 8169) 


2723 
(2855 l24l) 


STS(9) 
(9 x 12) 


0, 853, 217, 6, 2, 1108, 75, 20, 586, 1, 5 
0, 1797, 97, 3, 4, 485, 33, 37, 246, 1, 5 


5 


18 


(271760, 108704) 
((371100,92777) El) 


13588 


S-STS(13) 
(12 x 20) 


10484, 0, 12275, 10611, 9703, 10786, 10227, 11122, 3263, 7933, . . . 

3129, 21, 9554, 1, 2, 12183, 7837, 3084, 8297 
0, 0, 0, 12012, 13041, 498, 12534, 223, 7947, 8356, . . . 

12213, 13, 10701, 2, 1, 0, 9550, 10698, 4 



and continue with Step 2 until all c columns are used, 
that is, t = 0. 

By removing the last row and last K columns of the (J, ir- 
regular (c-6)xc base matrix B constructed using STS(c — b), 
we obtain a shortened (c — 6 — 1) x (c — K) (J, K — l)-regular 
base matrix B', which we denote -Bs-sts(c-&)- By deleting 
columns and rows, it is also possible to obtain intermediate 
codes, which are, however, irregular. 

Example 3: In the following we shall construct the (J = 3, K) 
base matrices B of dimension 9 x 12 (K — 4), dimension 
13 x 26 (A" = 6) and dimension 25 x 100 (K = 7). 
Using Algorithm STS, we obtain the following Steiner Triple 
Systems of order 9 (STS(9)), 13 (STS(13)) and 25 (STS(25)). 

STS(9) = 

{{2,3,5}, {1,4,6}, {1,3,7}, {2,6,7}, 

{4,5,7}, {1,2,8}, {5,6,8}, {3,4,8}, 

{1,5,9}, {2,4,9}, {3,6,9}, {7,8,9}} 

STS(13) = 

{{0,3,6}, {0,2,7}, {1,5,7}, {3,4,7}, 
{3,5,8}, {1,4,8}, {2,6,8}, {2,4,9}, 
{5,6,9}, {0,1,9}, {1,3, 10}, {0,4, 10}, 
{6, 7, 10}, {2, 5, 10}, {8, 9, 10}, {7, 8, 11}, 
{4, 6,11}, {1,2,11}, {0,5,11}, {3,9,11}, 



{10, 11, 12}, {7,9, 12}, {0,8,12}, {1,6,12}, 
{4,5,12}, {2,3,12}} 



{ 



STS(25) = 

{4,5,10}, 
{2,3,12}, 
{6,7,13}, 
{1,12,14} 
{5,6,15}, 
{0,8,15}, 
{2,8,16}, 



{1,9,10}, 
{0,9,12}, 
{0,5,13}, 
{0,2,14}, 
{3,10,15}, 



{7,8,11}, 
{6,8,12}, 



{1,6,11}, 
{8,9,13}, 



{2,10,13}, {3,4,14}, 
{7,9,14}, {5,11,14}, 
{4,12,15}, {1,7,15}, 
{6,10,16}, 



{11,13, 16}, {5, 7, 16}, 

{3,9,16}, {0,4,16}, {9,11,17}, 

{12, 13, 17}, {1,3, 17}, {4,7,17}, {0,6,17}, 

{2,5,17}, {8,17,18}, {3,11,18}, {2,4,18}, 

{13, 15, 18}, {0,10,18}, {1,16,18}, {6,14,18}, 



{9,18,19}, {4,8,19}, 
{0,3,19}, {2,7,19}, 
{17, 19, 20}, {9, 15, 20} 
{5,8,20}, {1,4,20}, 
{2,6,20}, {5,18,21}, {4,6,21}, 
{16,17,21}, {10,14,21}, {2,9,21}, 



{14, 15, 19}, {10, 11, 19}, 
{12,16, 19},{1,5, 19}, 
{10, 12, 20}, {0,11, 20}, 
{13, 14, 20}, {3, 7, 20}, 
{1,13,21}, 
{3,8,21}, 



{11, 15,21}, {7, 12,21}, {19, 21, 22}, {18, 20, 22}, 
{0,7,22}, {10,17,22},{3,5,22}, {6,9,22}, 
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{2,15,22}, {1,8,22}, {11, 12, 22}, {4, 13, 22}, 
{14, 16, 22}, {20, 21, 23}, {0, 1, 23}, {6, 19, 23}, 
{15, 16, 23}, {2, 11, 23}, {7,18,23}, {5,12,23}, 
{14,17,23}, {4,9,23}, {8,10,23}, {3,13,23}, 
{0,21,24}, {22, 23, 24}, {1,2, 24}, {16,20,24}, 
{7,10,24}, {8,14,24}, {13, 19, 24}, {3, 6, 24}, 
{12, 18, 24}, {15, 17, 24}, {5, 9, 24}, {4, 11, 24}} 

Each number 1,2,..., J occurs K times within the set of 
Steiner triples. However, the chosen Steiner triples are not 
uniquely determined. 

The corresponding base matrices of dimension 9 x 12 
STS(9), dimension 13 x 26 STS(13), and dimension 25 x 100 
STS(25) are sparse matrices with nonzero elements only in 
column i and row j, where the ith Steiner Triple contains the 
value j. The 9x12 (3, 4)-regular base matrix constructed from 
STS(9) denoted by -Bsts(9) is given, for example, by 



TABLE VII 

Properties of QC LDPC codes with girth g > 20 



B 



STS(9) 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 1 


12 


1 


fo 


1 


1 








1 








1 








o \ 


2 


1 








1 





1 











1 








3 


1 





1 














1 








1 





4 





1 








1 








1 





1 








5 


1 











1 





1 





1 











6 





1 





1 








1 











1 





7 








1 


1 


1 




















1 


8 

















1 


1 


1 











1 


9 


v> 























1 


1 


1 


1 / 



(23) 



Entries that correspond to edges in the base graph that can 
be, according to Algorithm STS, labeled with zero voltage are 
marked in bold. 

By removing the last row and the last K = 4 columns, 
the corresponding shortened 8x8 (3, 3)-regular base matrix 
_Bs-sts(9) follows directly as 



B 



S-STS(9) 



1 2 3 4 5 6 7 8 

/o 1 1 1 \ 

10 10 10 

1 1 

1 1 1 
1 1 
10 10 10 



1 



1 






111 



8 Vo 





111/ 



(24) 



This corresponds to removing the four Steiner Triples of 
STS (9) containing the number of the last row. Shortening the 
9 x 12 base matrix -Bsts(9) constructed from STS (9) to obtain 
a shortened 8x8 base matrix -Bs-sts(9) is unpractical as its 
code rate is R = 1 — 8/8 = 0. However, by shortening the 
13 x 25 base matrix Bsts(13) i n the same way we obtain a 
12 x 20 base matrix i?s-STS(i3) with the feasible code rate 
R = 8/20. 

In Table ED the obtained QC ( J = 3, iQ-regular LDPC 
block codes with girth g = 14, 16, and 18 constructed from 



K 


9 


(n, k) 


M 


Base graph (Table ITTA 


4 


20 


(1296000, 324002) 


36000 


(27 x 36), 3 = 8 


5 


20 


(31200000, 12480002) 


480000 


(39 x 65), 3 = 8 


6 


20 


(518400000, 259200002) 


4800000 


(54 x 108), 9 = 8 


4 


22 


(7200000, 1800002) 


200000 


(27 x 36), g = 8 (24) 


5 


22 


(325000000, 130000002) 


5000000 


(39 x 65), 3 = 8 


4 


24 


(39600000, 9900002) 


1100000 


(27 x 36), 3 = 8 



Steiner Triple Systems are presented. While the number of 
nonzero elements in each column is fixed to J = 3, the number 
of nonzero elements in each row K is specified in the first 
column. The second column corresponds to the obtained girth 
g, while in the third and forth columns we give the dimensions 
of the block code (n, k) after tailbiting to length M. And the 
fifth column contains which Steiner Triple System (STS(n)) 
is used. 

Finally, in the last column W we give the degrees of the 
corresponding degree matrix W in a compact way. As we 
have constructed the base matrices in such a way that the last 
nonzero entry in each column and the first entry in all other 
rows of the base matrix is labeled with a zero voltage, these 
entries are omitted. An entry of W' in column j and row i 
corresponds to the voltage degree of the (j + l)th nonzero 
entry in the ith row of the corresponding base matrix. 

Case 2-11: (J, K)-regular LDPC block codes 

When searching for QC (J = 3, if) -regular LDPC block 
codes with girth g = 20, 22 and 24, we started with previously 
obtained QC ( J = 3, if )-regular LDPC block codes of smaller 
block size and smaller girth, and (re-)applied our algorithms. 

The obtained results for QC ( J = 3, K )-regular LDPC 
block codes with girth g = 20, 22 and 24 are presented in 
Table IVHI They are all but one based on previously obtained 
( J = 3, K )-regular LDPC block codes with girth g = 8 
(cf. Table [TTTb. as specified in their last column in Table IVIII 
As before, the first column K denotes the number on nonzero 
elements in each column; then we give the obtained girth g 
and the dimensions of the block code (n, k) after tailbiting 
to length M. The corresponding degree matrices are too large 
and are omitted in Table IVIII but are available at 11391 . 

These codes are (probably) unpractical due to their huge 
block length. However, the table illustrates that by interpreting 
QC ( J, K )-regular LDPC block codes as base matrices and re- 
applying our algorithms we can find QC ( J, K )-regular LDPC 
block codes of "any" girth g. 

IX. Conclusions 

Using the relation between the parity-check matrix of QC 
LDPC block codes and the biadjacency matrix of bipartite 
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graphs, new searching techniques have been presented. Start- 
ing from a base graph, a set of edge voltages is used to 
construct the corresponding voltage graph with a given girth. 

By representing bipartite graphs in different ways, lower 
and upper bounds on the girth as well as on the minimum 
distance of the corresponding tailbiting block code have been 
discussed. 

New algorithms for searching iteratively for bipartite graphs 
with large girth and for determining the minimum distance of 
the corresponding QC LDPC block code have been presented. 
Depending on the given girth, the search algorithms are either 
based on all-ones matrices, Steiner Triple Systems, or QC 
block codes. Amongst others, new QC regular LDPC block 
codes with girth between 10 and 24 have been presented 
including their minimum distance if possible. In particular, 
the previously unknown minimum distance, for some known 
codes with girth 6 and 8, has been determined. 
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